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We study the onset and development of ledge instabilities 
during growth of vicinal metal surfaces using kinetic Monte 
Carlo simulations. We observe the formation of periodic pat- 
terns at [110] close packed step edges on surfaces vicinal to 
fcc(OOl) under realistic molecular beam epitaxy conditions. 
The corresponding wavelength and its temperature depen- 
dence are studied by monitoring the autocorrelation function 
for step edge position. Simulations suggest that the ledge 
instability on fcc(l,l,m) vicinal surfaces is controlled by the 
strong kink Ehrlich-Schwoebel barrier, with the wavelength 
determined by dimer nucleation at the step edge. Our results 
are in agreement with recent continuum theoretical predic- 
tions, and experiments on Cu(l,l,17) vicinal surfaces. 
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In surface growth under Molecular Beam Epitaxy 
(MBE) conditions, one of the central processes that con- 
trols the morphology of the surface is mass transport be- 
tween growing layers. In homoepitaxial growth, there is 
usually an additional energy barrier controlling the inter- 
layer transitions, known as the Ehrlich-Schwoebel (ES) 
or step edge barrier []l| . It plays a central role in stabiliza- 
tion of certain growing crystalline facets, and eventually 
leads to growth of mound-like structures with a dynami- 
cally selected slope and length scale . The influence of 
an interlayer ES barrier on growth has been extensively 
studied and its role on growth is now well understood 

However, recently it has been realized that in 1 + 1 
dimensional ledge growth corresponding to step-flow ge- 
ometry, there is an analogous phenomenon which is due 
to the additional kink Ehrlich-Schwoebel energy barrier 
for going around a kink site at the step edge. The cor- 
responding kink Ehrlich-Schwoebel effect (KESE) gen- 
erated by it leads to the growth of a regular instabil- 
ity at step edge with a dynamically selected wavelength 
[g[. This is in contrast to the so-called Bales-Zangwill 
instability (BZI) |^ which tends to destabilize the ID 
ledge morphology during growth because of terrace dif- 
fusion and step crossing, with no assumptions about line 
diffusion along the ledge. Ledge instabilities were orig- 
inally found and reported experimentally on Cu(l,l,17) 
vicinal surfaces 0[ but attributed to the BZI scenario. 
More recent experiments on the Cu(l,l,17) surface pro- 
pose that the KESE instability may lead to formation 



of regularly shaped patterns with dynamical wavelength 
selection Recent theoretical studies of such instabili- 
ties suggest that KESE may indeed supersede BZI in the 
formation of growth patterns [^,^ . 

The growth of instabilities and wavelength selection 
have been studied recently within the framework of a con- 
tinuum step model [^ and a simple solid-on-solid (SOS) 
lattice model For the latter case, it has been demon- 
strated by kinetic Monte Carlo (KMC) simulations that 
KESE leads to the formation of wavy steps [^. However, 
so far there is neither detailed knowledge of the actual 
structure of the patterns nor their dynamical evolution 
under realistic MBE conditions. The predictions of the 
continuum theory are mainly given in the weak KESE 
limit but they are not directly applicable in interpreting 
the experiments [^,|l6|. Namely, for simple metals such 
as Cu surfaces vicinal to fcc(OOl), at close packed step 
edges a strong KESE is expected. This is because theo- 
retical estimates indicate barriers of the order of 0.5 eV 
for jumps around the kink site in the close packed [110] 
direction [|l^. This is nearly twice the barrier of 0.26 eV 
for jumps along the straight step edge [|l^. Based on the 
symmetry between the interlayer ES barrier which van- 
ishes for [100] step edges [0, it is expected that the kink 
ES barrier vanishes for these orientations, too. Therefore, 
on surfaces vicinal to fcc(OOl) one should see a clear dif- 
ference in unstable growth between the [100] steps and 
the [110] steps, the latter case being dominated by KESE. 

In this letter we present a detailed study of the growth 
and morphology of the steps on the Cu(l,l,17) surface 
vicinal to Cu(OOl). The model system used here is based 
on KMC simulations of a lattice-gas model [^ with ener- 
getics obtained from the effective medium theory (EMT) 
potential As discussed in detail in Ref. [|o| the EMT 
barriers and their relative ordering is in good agreement 
with available experimental data for the Cu(OOl) surface. 
The intralayer hopping rate v of an atom to a vacant 
nearest neighbor (NN) site can be well approximated by 
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vo exp [-l3{Es + min(0, 1^nn)Eb] 
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Fig. 2: Snapshots of typical step edge profiles with 

step orientations in the (a) and (b) [flO] and (c) [100] 
directions at T = 240 K with F = 6 x 10"^ ml/s fo^. 
coverages 6 —1.0, 5.0, and 4.0, respectively. In the latter 
case KESE vanishes and the instability is strongly sup- 
pressed. The lateral and vertical scales are given in units 
of the lattice constant. 



Fig. 1: Ledge geometry and some relevant single atom 
jump processes on (a) the close packed [110] step edge 
and on (b) the open [100] step edge. The energy barriers 
are given in eV. 

where the attempt frequency fo = 3.06 x 10^^ s~^ and 
the barrier for the jump of single adatom is Es = 0.399 
eV. When there is at least one atom diagonally next to 
the saddle point the barrier is Es = 0.258 eV. The change 
in bond number — 3 < Ann < 3 is the difference on the 
number of NN bonds between initial and final states with 
bond energy Eb — —0.260 eV. The interlayer processes 
are also included in the model. Step crossing has barrier 
of 0.58 eV from a straight edge and 0.44 eV through the 
kink site ||l^. We note that within the EMT, barriers 
on the Ag(OOl) and Ni(OOl) surfaces are very similar to 
those of Cu(OOl) up to an overall scaling factor ||l^. The 
geometry of the step edges with some of the pertinent 
single jump processes are shown in Fig. 1. The energetics 
of the model give a kink barrier of 0.518 eV for the [110] 
steps, and vanishingly weak KESE with a barrier of only 
0.002 eV for the open [100] steps. 

The KMC simulations were implemented using the al- 
gorithm by Bortz et al. |13| with a binary tree struc- 
ture [0. This allowed us to grow up to five mono- 
layers of Cu under realistic temperature and flux con- 
ditions to study the development of unstable ledge pat- 
terns. The explored temperature range was T — 240—310 
K and the flux F = 3 x 10"^ _i q mL/s. Thus the ra- 
tio between the terrace diffusion and the deposition flux 
D/F ^ 6 X 10^ — 9 X 10'' in units of the lattice constant 



a = 0.255 nm, corresponding to a typical MBE regime 
p^ . System sizes considered consisted of eight terraces, 
each of width Lj. — 8.5 and length up to Ly — 2000. 
With fully periodic boundary conditions there was a de- 
pendence of the selected wavelength on Ly. With open 
boundaries this dependence disappears, and in the ledge 
direction there was no difference in the results between 
Ly — 500 and Ly = 1000. However, periodic boundary 
conditions were necessary to study the temperature de- 
pendence of the structures, in which case we explicitly 
checked that this feature did not depend on Ly. 

The choice of the initial conditions warrants some dis- 
cussion. In principle, ID step edges are thermally rough 
in equilibrium. Thus, we chose to start from initially 
rough step edges. We also performed some test runs 
starting from ideal, smooth ledges. In these cases devel- 
opment of the instability was much slower than when the 
steps were rough. We also discovered slight dependence 
of the characteristic wavelength of the instability on the 
initial condition. For both cases, however, we verified 
that the same temperature dependence was obtained. 

We first show results for the step orientation [110] 
where KESE is relatively strong. In Figs. 2 (a) and 
(b) there are step profiles after 6* = 1.0 and 5.0 mono- 
layers, respectively, at T = 240 K with F = 6 x 10~^ 
ML/s. The development of an instability with a rela- 
tively well-defined wavelength for all ledges is apparent, 
and the steps seem to "phase-lock" at the largest cov- 
erages studied. To quantify these results, we consider 
the lateral height correlation function for step edge pro- 
files C,{x,9) = h{x,9) — h{9), where h is the step height 
measured from y = 0, and h is its spatial average: 
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Fig 3: The correlation function C{x, 9) for step edge 
fluctuations along the [110] ledge, for coverages 9 =0.4, 
1.0 and 2.0 at T = 270 K. The corresponding Fourier 
components of C are shown in the inset. An estimate 
for the wavelength obtained from the average distances 
between fingers is denoted by an arrow. 

Here, the brackets denote an average over the (deposi- 
tion) noise and over all ledges in the system. The corre- 
lation hmction C{x^ 9) for the [110] step edges at T = 270 
K with F = 6 X ML/s is shown in Fig. 3 for 

9 — 0.4, 1.0, and 2.0, and its Fourier components are 
shown in the inset. The first minimum of C{x, 9) gives a 
measure of the characteristic length scale A « 37 nm in 
the system, and as can be seen from Fig. 3 it is almost 
independent of 9 even for coverages below the apparent 
phase-locking. The same information is contained in the 
Fourier components of ^. The arrow in Fig. 3 shows 
an estimate of the average distance between the growing 
fingers, as obtained directly from the ledge configura- 
tions. We note that in experiments the stabilization and 
phase correlation of the perturbations seems to improve 
between 5 ML to 50 ML [^. Due to computational re- 
strictions we have not been able to probe this regime, 
however. 

In Fig. 2 (c) we show some step profiles for the open 
[100] steps after 9 = 4.0 monolayers. The difference 
with respect to the close-packed steps is striking. Within 
the accuracy of the data, we find no evidence of unsta- 
ble growth but the ledge fluctuations appear completely 
random. This is strong evidence indicating that in the 
present case, the patterns on the [110] ledges are due to 
KESE and not by BZI. In the absence of KESE, enhanced 
edge diffusion along the sides of the steps here strongly 
suppresses fluctuations of the [100] ledges (cf. Fig. 1 
(b)). 

Let us next consider the origin of the wavelength 
(A w 30 nm at T = 240 K with = 6 x 10"^ML/s) ob- 
served for the [110] ledges. It is of the same order of mag- 
nitude as the wavelength obtained in the experiments [pj . 
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Fig. 4: Temperature dependence of the the wavelength 
A obtained from the distance between the fingers. It has 
been determined after deposition of 2.0 ML with F = 
6 X 10-3 ML/s at T = 240 - 310 K. The least squares 
fit is shown by a continuous line with a slope denoted in 
the figure. In the inset the flux dependence of A is shown 
at T = 300 K with fluxes F = 3 x 10"=^ - 1.0 ML/s after 
deposition of 2.0 ML. The straight line is the theoretical 
prediction of Eq. (3). For the highest flux values, there 
is indication of terrace nucleation. 

According to continuum theory for KESE [|^, there are 
two important length scales controlling step flow growth 
with KESE. They are denoted by for dimer nucleation 
on a straight step, and the the kink Schwoebel length 

= exp[{Es - Ed)/kBT] - 1 §], which is related to the 
energy barriers Es and Ed for jumps around a kink site 
and along a straight edge, respectively. When is ^ ^c, 
the nucleation length can be estimated from the model 
of Politi and Villain jl^ (see also Ref. [^). However, for 
the [110] ledges, ^ 10^ and 4 ^ 10^ around room 
temperature, which indicates that we must consider the 
opposite case of a strong KESE. 

In the case of a strong KESE, a more appropriate 
estimate for the nucleation length £c is obtained by 
considering the probability for nucleation of dimers at 
straight step edges. In the limit of infinite KESE it 
has been shown that 4 ~ (Ds/Fs)'^^^ [0- Here Ds 
oc exp[— Fd/fesT] is the diffusion coefficient along the 
straight ledge and Fs = FLx is the fiux at the ledge. In 
the case of finite but large £s, the corresponding ID diffu- 
sion equation with appropriate boundary conditions has 
been solved by Politi jlj]. From his solution, assuming 
that £c is controlled by dimer nucleation at the bottom 
of the growing step structure gives ic = (12Ds/Fs)-^/*. 
Since the width of the ID terrace fluctuates between zero 
and £c J0, we obtain an estimate for the wavelength as 
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The scaling exponent 1/4 is a direct consequence of dimer 
nucleation at the ledge and is always expected for strong 
KESE within the present model. The relation in Eq. (3) 
predicts that the temperature dependence of A is con- 
trolled by an effective barrier E^s = Ed/ 4:. In the present 
case this becomes about 65 meV. The estimate for Ecs 
in the case of BZI gives a value nearly an order of magni- 
tude larger Eq. (3) also predicts that A ~ with 
a — 1/4, in contrast to the BZI case where a = 1/2 

With the present energetics, Eq. (3) gives A « 40 nm 
at T = 270 K with F = 6 x 10"^ ML/s which is in very 
good agreement with the simulation result A « 37 nm. 
This is in contrast to BZI which yields A « 1 nm We 
have furthermore tested the predictions of Eq. (3) by es- 
timating A at different temperatures and also by varying 
the flux. In Fig. 4 we show the temperature dependence 
of A for various temperatures in an Arrhenius plot. A 
straight line fit to the data gives Ecff — 75 ± 10 meV, 
in good agreement with the prediction based on Eq. (3). 
This estimate is also in good agreement with the experi- 
mental result by Maroutian et al. In the inset of Fig. 
4 we also show the flux dependence of A and compare 
with that obtained form Eq. (3) with no fitting param- 
eters. We find excellent agreement between theory and 
simulations, except for the largest fluxes where terrace 
nucleation becomes apparent. Fitting to flux values up 
to F 9 X 10^2 ML/s giygg Q, = 0.23 ± 0.03. 

Finally, we discuss the roughness of the step in terms of 
the width of the step edge fluctuations 'w{9) = \/ {({Oy). 
We find that w{d) increases with coverage and does not 
show signs of saturation up to about 5 ML. Up to this 
region the width follows a power law behavior w{6) ~ 9^ , 
with f3 « 0.3 as seen in case of isolated step [|l9|. The 
solid-on-solid model gives for strong KESE an exponent 
/3 « 0.57 [D, and /3 = 1 has been obtained for a one-sided 
step growth model with infinite kink barrier This 
indicates that f3 depends on the details of the model for 
unstable growth. 

To summarize, the general situation regarding ID 
ledge instabilities under growth is a complicated one ^] 
(see also |^,|o|). The simulations presented here give 
support to the view that on vicinal Cu(l,l,m) surfaces 
the observed instability is due to KESE. This means that 
the competing BZI is of no importance neither in length 
and time scales studied here nor in those accessed in the 
experiments. In our simulations we have been able to 
demonstrate the onset of the wavelength selection, in 
qualitative agreement with continuum theoretical predic- 
tions. Moreover, it was possible to connect the observed 
temperature dependence of wavelength to the underly- 
ing energetics. This is in qualitative agreement with ex- 
perimental observations and theoretical predictions for 
strong KESE, and clearly incompatible with BZI. Simu- 
lations also confirm the crucial role of the dimer nucle- 
ation length in determining the length scale of the pat- 
terns. Although we cannot conclusively demonstrate the 



phase-locking of the step fluctuations, there is evidence 
that this indeed occurs at higher coverages. 
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